function SAT_SST_Step_02_calculate_anomalies_for_raw_SAT(P)

    % *********************************************************************
    % Load Scaled Air tempertures
    % *********************************************************************
    file_load = [SAT_SST_func_IO('data',P),'/Raw_SAT_SST_',P.date,'.mat'];
    load(file_load,'raw_SAT_full','lon','lat');

    % *********************************************************************
    % Call external function to connect stations and calculate anomalies
    % *********************************************************************
    % use the most direct way to compute climatology, no fitting!
    [raw_SAT_anm,Tier] = SAT_SST_func_connect_stations(raw_SAT_full,lon,lat,...
                            P.do_gradual_match,P.distance,P.cal_anomalies);

    % *********************************************************************
    % Save Data
    % *********************************************************************
    file_save = [SAT_SST_func_IO('data',P),'/Raw_SAT_anm_',P.date,'.mat'];
    save(file_save,'raw_SAT_anm','Tier','-v7.3');
end
